TITLE 

TYPE 

METHOD OF USE 



CAPACITY 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



LIBRARY ROUTIN E K8- 189 
By Gene Ho Golub 



Product Moment Correlations, Variance-Covariances, 

Means and Standard Deviations for Use with Magnetic Drum 

(SADOI Only) 
Entire program 

1) Read in the master program. A sum check is performed. 
If the program is read correctly, the machine will stop 
on jk OIK; otherwise on an FF 01? order. 

2) Read in the parameter tape. The computer stops on 2k QJ$, 

3) Read in the data tape. 

k) After all data have been read into the computer, the 
desired information will be printed out. The program stops 
on 2k OIK. 

5) Another problem can be begun by reading in new parameters. 
The maximum number of variables is lkk$ there is no limit on 
the number of observations. 



PUNCHING OF THE TAPE 



For every problem four parameters are necessary. They are 
as follQws and may be punched in any order on the tapes 
l) L^i"f be the number of decimal places to which the 
cot*re|aiiQii inatrix is to be printed. Put fF on the parameter 
i&pe. If no print out is desired, f = 0. 
2 J Let "a- m be the number of decimal places to which the 
Variance covariance matrix is to be printed Put JlL oh tjie 
parameter tape. If no print out %s desired, JC= 0. 
3) Let "j" be the number of decimal places to which the 
means and standard deviations are to be printed. |*ut JJ 
pj* the. parameter tape. If no printout is desired, j = 0. 
k) Let "s* be the size Of the sample. Put sS on the 
parameter tape. 



- 2 - 

This must be followed by either a "0" or "1". A "0" 
indicates that only a triangular matrix is to be punched 
row by row, while a 1 indicates that the square matrix 
is to be punched. 

Each observation (which must lie in the range - 1 < x < l) 
is punched as a sign followed by up to 12 decimal digits. 
The character N must be punched after each row. If an 
F follows a row, then the machine will stop and another 
part of the data tape may be placed in the reader. 

THE ERJHT OUT The correlation matrix, scaled by one-tenth, is printed 

in the first column. The covariances appear in the 
second column. If the full matrix is printed out, then 
an N appears after each column of the, matrix. Following 
this the means and standard deviations are printed in two 
parallel columns. 

NOTES i) Correlations with constants are assumed to be zero. 

In order to avoid a division hangup the correlation 
between a constant and itself will be zero. 

2) If the incorrect number of variables has been punched 
in a row, then the machine will hangup on an FF 015 order 
at location (O^N)^. 

3) If the variance is negative and outside the tolerance 
limits, then the machine will hangup on an FF 0l6 order at 
(088).. c, otherwise it is set equal to zero. 

k) An FF023 stop in location (081)^ indicates an 
arithmetic error, possible from scaling. 



DAT E March 19, 1957 BT: H/27/59 

PROGRAMMED B Y G. H. Golub 

APPROVED BY D. Eo Muller ______ 



ns 



TIME ESTIMATES FOR K-8 PROGRAM 

The following times are estimates only They are for use as a guide to help provide closer estimates in 
the time required to run K-8 programs through Illiac. In general, estimates based on these times should not 
be more than + 10$ in error. 

The chart below gives the time for fixed intervals of variables in seconds per sample, i.e., for a 100 
sample size problem, multiply the time in column two by 100 to obtain the time in seconds the problem is to 
run before output occurs. The remainder of the chart gives the time, in minutes , that a square or diagonal 
matrix would require to output. For convenience sake, these times are given for 3/ 4, 5> and 6, decimal places 
of output. 







READ-IN TIME 


PRINT-OUT TIME IN MINUTES 










NO. 


OF VARIABLES 


SEC/SAMPLE 


TRIANGULAR MATRIX 




SQUARE MATRIX 












3 dec. 4 dec. 5 dec 6 dec. 


3 dec. 


4 dec. 5 dec. 


6 dec. 






5 


1 


1111 


1 


1 1 


1 






10 


1 


1111 


1 


1 1 


1 






15 


2 


1111 


1 


1 1 


1 






20 


2 


111 1 


1 


1 1 


2 






25 


3 


1111 


1 


2 2 


2 






# 


4 


1 1.5 1.5 2 


2 


3 3 


4 






35 


5 


1.5 2 2 2.5 


3 


4 4 


5 






40 


6 


2 2.5 2.5 3 


4 


5 5 


6 






45 


7 


2.5 3 3 3.5 


5 


6 6 


7 






50 


7-5 


3 4 4 4.5 


6 


8 8 


9 






55 


9 


4 4.5 5 5.5 


8 


9 10 


11 






6o 


11 


4.5 5.5 6 6.5 


9 


11 12 


13 






65 


13 


5.5 6.5 7 7-5 


11 


13 14 


15 






70 


14 


6.5 7.5 8 9 


13 


15 16 


18 






75 


16 


7.5 8.5 9.5 10.5 


15 


17 19 


21 






- ab" " 


ia 


8.5 9.5 10o5 11.5 


17 


19 21 


23 






85 


20 


9.5 11 H.5 13 


19 


22 23 


26 






90 


22 


10.5 12 13 14.5 


21 


24 25 


29 






95 


24c5 


11.5 13.5 15 16 


23 


26 29 


32 






100 


27 


13 15 17 18 


25 


29 33 


35 






105 


31 


14 16.5 18 19.5 


27 


32 35 


38 






110 


34 


15.5 18 20 21.5 


30 


55 39 


42 






115 


37 


17 19.5 21.5 23o5 


33 


37 42 


46 






120 


41 


18.5 21 23.5 25.5 


36 


41 46 


50 






125 


44 


19.5 23 25o5 27.5 


38 


45 50 


54 






130 


48 


21 25 27o5 30.5 


41 


49 54 


60 






1^5 


52 


22.5 26o5 29.5 33 


44 


52 58 


65 






140 


56 


25 28.5 31.5 35 


>9 


56 62 


69 






144 


61 


27o5 31 34 38 


5 1 * 


61 67 


75 





(M. Foster) 



LOCATION 



3 

k 

5 
6 

7 
8 

9 
10 



ORDER 



NOTES 



PAGE 1 



MA 



00 3K 
00 F 
00 (S3) 

00 F 
00 217F 
00 F 

00 2560F 

00 F 
00 (F2) 
00 F 
00 (Rl) 

00 12K 

50 OF 

kO 102^F 

06 11F 

00 S5 

00 IF 

00 IF 

50 OF 

L5 1024F 

05 11F 

00 S5 

00 F 

00 IF 

00 F 

00 1000 0000 0000 J 

00 F I 

00 10F I 

kO (S3) I 
i 

00 (S3) I 
80 F 
00 F 
00 F 
00 F 



LOCATION 
11 



12 



13 



3 

k 

5 
6 

7 
8 

9 

10 

11 

12 

13 



ORDER 
00 F 

00 F 
JO F 
kO F 
JO F 
L5 F 
00 26K 
19 3F 
kO 2F 
92 151F 
kl F 
81 4F 
LO 19F 
32 6L 
IA 19F 
50 OF 
7^ 19F 
S5 F 
kO F 
26 2L 
42 7L 
L5 OF 
kO F 
L5 2F 
IA 2F 
kO 2F 
32 1L 
IA IF 
40 9F 
81 kF 
00 39F 
40 10F 
19 38F 
66 IF 
S3 F 



NOTES 



PAfig g 



r 

K 8 



from 9 
from 6 



Read in parameters and convert 



from 3 



by 6 



LOCATION 
14 

15 
16 

17 
18 

19 



ORDER 



21 
22 

23 
24 

25 
26 

27 
28 
29 
30 
31 



1+0 

4l 
50 
50 
24 
LO 
30 
L5 
40 
LD 
46 
10 
42 
42 
40 

L5 
L4 
40 
L4 
40 
L4 
42 
42 
L5 
50 
74 

S5 
10 

L4 

L4 
40 

L3 
40 

L5 
42 
41 



7F 

6f 

(S3) 

15L 

(N12) 

(2) 

17L 

21S4 

11F 

20F 

6F 

20F 

6F 

21F 

OF 

20F 

6F 

22F 

6f 

23F 

6f 

24F 

25F 

OF 

OF 

OF 

F 

IF 

OF 

OF 

OF 

24F 

32L 

22F 

32L 

IF 



NOTES 



PAGE 3 K 8 



from 17 



l/s 



Read in first row 



Count hov many variables 



mm 



LOCATION 



ORDER 



NOTES 



PAGE k 



K 8 



32 


JO OF 




kO F 


33 


F5 32L 




40 32L 


3^ 


LO 12F 




36 36L 


35 


l4 13F 




kO 32L 


36 


F5 IF 




40 IF 


37 


LO OF 




36 39L 


38 


27 32L 




23 161L 


39 


L5 kF 




00 20F 


ko 


46 161L 




L5 5F 


ki 


00 20F 




k6 l6kL 


k2 


L5 3F 




00 20F 


^3 


k6 183L 




46 185L 


kk 


ki 8f 




22 50L 


k5 


50 (S3) 




50 1+5L 


k6 


26 S4 




LO (2) 


kj 


30 47L 




L5 21S4 


48 


LO 11F 




40 OF 


k9 


L3 OF 




32 50L 



from 38 
hy 30,31,33,35 



Clear enough store for data 



from 34 



from 37 



from 97 



from 47 



Read in data 



Are number variables the same? 



LOCATION 




52 
53 
5 4 

55 



ORDER 



FF 21F 
L5 22F 
42 55L 
46 55L 
L5 20F 
42 54L 
46 T6L 
4l 2F 
50 7F 
7 J F 
L4 F 
40 F 



F5 54L 





42 54L 


57 


L5 55L 




L4 14F 


58 

1 


40 55L 




F5 2F 


59 


40 2F 




LO 21F 


60 


56 54L 




40 OF 


61 


L5 25F 




40 69L 


62 


L5 24F 




40 88L 


63 


4l IF 




L5 20F 


64 


42 76L 




'L5 23F 


65 


42 70L 




42 79L 


66 


46 79L 




| 42 87L 


67 


1 F4 OF 




\ 42 133L 



NOTES 



PAGE 5 



K8 



from 44, 49 



from 60 
by 52, 56 
by 51, 58 
by 51, 58 



J=l 



X ij /S 



P = 1, 2, 
i = 1, '2, 



n 



from 95 



S LOCATION 


ORDER 


NOTES PAGE 6 K 8 


1 68 


F5 OF 




Read row of matrix into working space 




40 OF 






69 


JO OF 


from 75 






L5 F 


by 61, 71, 73 




TO 


JO OF 








40 F 


by 65, 74, 130 




71 


F5 69L 
1+0 69L 






72 


LO 15F 








36 74L 






73 


L4 i6f 








ko 69L 






!**■ 


F5 70L 


from 72 






42 70L 






75 


LO 133L 








36 69L 






76 


50 P 


by 53,855 from 84 






7J F 


by 64, 80 




77 


kO 2F 

50 2F 




x ik x jk/ s $ = *»' 2 * ••"' x 
i, j = 1, 2, ..., n 


78 


7J 7F 

36 79L 






79 1 


L4 F 


by 66, 82 






40 F 


by 65, 82 




80 


F3 76L 








k2 76L 






8l 


L5 79L 
L4 l4F 






82 


kO 79L 
F5 IF 






83 


kO IF 
LO OF 






84 


36 76L 
L5 76L 






85 


L4 l4F 








k6 76L 







LOCATION 


1 

ORDER 


NOTES PAGE 7 K 8 

- 


86 


L5 87L 








L4 OF 






87 


40 2F 




■ 




L5 F 


| by 66, 92$from 93 


88 


JO OF 








40 F 


by 62, 89, 91 




89 


F5 88L 




' 




40 88L 




Store row back into 


90 


LO 12F 
36 92L 




matrix 


91 


L4 13F 








40 88L 






92 


F5 87L 


from 90 






42 87L 


! 




93 


LO 2F 


! 






32 87L 


i 




9 k 


L5 IF 
LO 21F 




' 


95 


36 63L J 


\ 






F5 8F 1 






96 


40 8F 








LO 9F 






97 


36 45L 


■ 






L5 25F 






98 


40 102L 
L5 20F 


i 
1 




99 


42 111L 


1 






L5 22F 






100 


42 105L 1 
46 105L 


. 




101 


4l OF 








4l 11F 






102 


JO OF 


from 119 


- 




L5 F 


by 98, 113, 115 i 








1— _ 



LOCATION 
103 

104 
105 
106 
107 
108 

109 
110 
111 
112 

113 
114 

115 

116 

117 
118 

119 
120 



ORDER 
36 104L 
FF 35F 
40 IF 
00 IF 
50 F 
79 F 
L4 IF 
40 IF 
L3 IF 
L4 9F 
36 109L 
23 111L 
L5 IF 
32 110L 
FF 22F 
50 110L 
22 (Rl) 
40 F 
F5 11F 
F4 102L 
40 102L 
LO 15F 
32 115L 
L4 loF 
40 102L 
F5 111L 
42 111L 
L5 105L 
L4 14f 
40 105L 
F5 11F 
40 11F 
LO 21F 
36 102L 
40 6F 
L5 25F 



NOTES 



PAGE 8 



K8 



from 103 



"by 100, 117; from 114 
by 100, 117 



from 108 



from 109 



by 99, 116; from 108 



Ex 2 



L/» - '4 



k=l 



Get standard deviation 



LOCATION 



ORDER 



NOTES 



PAGE 9 



K 8 



121 
122 

123 
12^ 

125 
126 

127 
128 
129 
150 
131 
132 
133 
13^ 
135 
136 

137 
138 



from 179 



from 138| 1U7 

by 127,13^136,1^1^6" 



by 67, 137 



from 135 



Read in row of matrix 
into working space 



I LOCATION 
j 139 

I i4o 
« 

l 

| l4l 
t 

l42 
143 
144 
I l45 
146 
147 
148 
149 
150 
151 
152 
153 
154 

155 
156 



ORDER 

L3 8F 
32 141L 

41 8F 
L5 132L 
40 JF 
L5 10F 
32 148L 
L5 132L 
L4 11F 
LO 17F 
40 132L 
LO 15F 
32 146L 
L4 16F 
40 132L 
F5 11F 

40 11F 
FO 21F 
36 132L 
L5 20F 

42 I56L 
L5 23F 
46 154L 

41 8F 
L5 22F 

42 15 3L 
92 131F 
92 515F 
50 F 

79 F 
L4 F 
40 9F 
L3 4f 
32 162L 
50 F 
7J F 



NOTES 



PAGE 10 K 8 



from 139 



from 145 



from 142 



from 170; 173 

by 123, 178 
by 151, 166 
by 149, I67 



by 122, 177 
by 148, 168 



E 

k=l 



x., x„ /s - X. X. 

ik jk' 1 3 



Compute product of standard deviation 



LOCATION 


ORDER 


NOTES PAGE 11 K 8 | 


157 


kO OF 








L3 OF 




t 


158 


32 38L 




1 




50 9F 






159 


75 18F 








66 OF 




■ 

! 
| 


160 


S5 OF 






1 


32 161L 






16X 


50 F 
50 161L 


TDy *K) 


Print out correlations 


162 


26 (P2) 




. 




L3 5F 


from 155 




165 


32 165L 








L5 9F 






164 


50 F 
50 l6^L 


Sy M 


Print out variance 


165 


26 (P2) 








F5 153L 


from 163 




166 


k2 153L 




. 




| L5 15^L 






167 


IM lh? 








k6 15^L 






168 


F5 156L 








te 1&&L 






169 


F5 8F 




; 




kO 8F 






170 


IO 6F 








36 152L 






171 


L5 lOF 
32 17^L 






172 


L5 8F 








K) 21F 






1T5 


36 152L 
92 1J1F 






17k 


92 TK® 








L5 6F 


trm 171 





; LOCATION I 


ORDER 




NOTES PAGE 12 K 8 


i f 

175 | 


LO 21F 




: 


t 


32 179L | 




1 


176 


L5 156L 








L4 l4F 






177 


46 156L 








L5 153L 






178 


L4 l4F 
46 153L 






179 


26 126L 








L3 3F 


from 1255 175 




180 


32 190L 
92 139F 






181 


4l 5F 








92 131F 




Print out means 


182 


92 515F 








L5 F 


by 124, 187 




183 


50 F 


by 43 






50 183L 






184 


26 (P2) 




Print out standard deviations 




L5 F 


by 122, 188 




185 


50 F 
50 185L 


by 43 




186 


26 (P2) 

F5 182L 






187 


42 182L 
F5 184L 






188 


42 184L 
F5 3F 






189 


ho 3F 








LO 21F 


■ 




190 


32 181L 








24 L 


from 180 




(N12) 


00 K 

00 K(P2) 








Routine P2 


- 52 


Print (A) with or without Sign to n 
Places as Determined by a Program 






Parameter 



LOCATION 

(2) 
(S3) 



ORDER 

00 K 
00 F 
00 2F 
00 F 
00 F 

00 800K 
80 40F 
LO F 
ifO F 

L3 F 
3k 26F 
FF 23F 
26 800N 



NOTES 



I 

I 78N04iK)N070 



PAGE 13 



K 8 



Sum Check 



